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We discuss the Hawking temperature of near-equilibrium black holes using a semiclassical analysis. 
We introduce a useful expansion method for slowly evolving spacetime, and evaluate the Bogoliubov 
coefficients using the saddle point approximation. For a spacetime whose evolution is sufficiently 
slow, such as a black hole with slowly changing mass, we find that the temperature is determined 
by the surface gravity of the past horizon. As an example of applications of these results, we study 
the Hawking temperature of black holes with null shell accretion in asymptotically flat space and 
. . . the AdS-Vaidya spacetime. We discuss implications of our results in the context of the AdS/CFT 

04 , correspondence. 
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I. INTRODUCTION 



The fact that black holes possess thermodynamic properties has been intriguing in gravitational and quantum 
theories, and is still attracting interest. Nowadays it is well-known that black holes will emit thermal radiation with 
Hawking temperature proportional to the surface gravity of the event horizon This result plays a significant 
role in black hole thermodynamics. Moreover, the AdS/CFT correspondence opened up new insights about 
(jj I thermodynamic properties of black holes. In this context we expect that thermodynamic properties of conformal field 
f~i . theory (CFT) matter on the boundary would respect those of black holes in the bulk. In particular, stationary black 
holes correspond to thermal equilibrium of the dual field theory at finite temperature equal to the bulk Hawking 
. temperature 0]. 

^ ■ It is interesting to study if such thermal properties persist when time dependence is turned on. From a technical 
\ viewpoint, the original derivation of the Hawking radiation [l[ was performed on a static or stationary background 
OO ■ (precisely speaking, as an approximation at late time) , and it is not obvious how to accommodate dynamical spacctimcs 
^^O [ into the scheme. From a physical viewpoint, such a study can be useful since realistic black holes in our universe are 
. dynamical or at least quasistationary. For example, processes like black hole formation after a gravitational collapse 
or black hole evaporation due to the Hawking radiation may involve highly time-dependent phases, and it is not 
obvious how the Hawking radiation will behave in such systems. 

When time dependence is sufficiently weak, however, we may expect thermodynamic properties of a black hole 
to persist. This expectation partially stems from physical insight on ordinary thermodynamics systems, in which 
thermodynamic properties persists when the system is sufficiently near equilibrium and quasistationary. If black 
hole thermodynamics is robust enough, such insight on usual thermodynamic systems leads to the above expectation 
^ ^ . on dynamical spacetimes and black holes. Many studies have been made on generalizing Hawking temperature for 
■ dynamical black hole approaches (3-[l0j. and they support such an expectation. 

Within the context of the AdS/CFT correspondence, thermodynamic properties of dynamical black holes lead to 
new insights into the dynamics of strongly coupled field theory on the boundary. One typical problem in this line 
is thermalization processes in the boundary field theory, which is holographically modeled by formation of the bulk 
black hole and its equilibration into a stationary state. Much effort is devoted to such problems in connection with, 
for example, application of the AdS/CFT correspondence to QCD physics like that in the RHIC experiment, or to 
nonequilibrium phenomena in condensed matter physics and fiuid dynamics pdl - |23j . In this context, the Hawking 
radiation in the bulk is interpreted as quantum fiuctuation in the boundary theory, and it plays an important role in 
certain setups [13, HH . 

Based on these interests, in this paper we consider nonstationary spacetimes, and try to determine the Hawking 
temperature measured by observers in the asymptotic region. One naive way to do this would be to define a certain time 
coordinate to associate temperature of the black hole (future) horizon to that measured by asymptotic observers. Even 
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though this procedure is seemingly natural from the viewpoint of black hole thermodynamics or horizon dynamics, 
it is puzzling from the viewpoint of causality, since this temperature is determined using information in the future 
region for that observer. In the AdS/CFT setup, if we were to define temperature at a point on the boundary in 
this way, we would need information in the bulk region, which is not causally accessible from that boundary point. 
Such a property would not be desirable if we care about causality in the bulk and boundary, especially when time 
dependence comes into the story. 

To avoid such problems, we consider the conventional derivation of the Hawking radiation based on Bogoliubov 
transformations for nonstationary background. An advantage of this method is that we can naturally determine 
temperature measured by asymptotic observers using only information causally available to them.^ Particularly, we 
focus our attention on an eternal black hole rather than black hole formation by gravitational collapse because we are 
interested in near-equilibrium system such as transition from an equilibrium state to another equilibrium state. ^ In 
other words, we focus on late-time dynamics after the black hole has formed. We clarify how the Hawking temperature 
changes when a black hole becomes dynamical, for instance, due to mass accretion to the black hole, and show that 
the Hawking temperature of dynamical black holes can be naturally associated with "surface gravity" of the past 
horizon. 

The paper is organized as follows. In Sec. [Hi we evaluate the Bogoliubov coefficients using the saddle point 
approximation and show that the "surface gravity" of the past horizon gives the Hawking temperature observed 
in the asymptotic region. In Sec. IIIIl as illustrations of the use of our method, we consider applications to simple 
examples in asymptotically flat and AdS cases. We summarize and discuss implications of the results in Sec. IIVI 



II. ESTIMATION OF BOGOLIUBOV COEFFICIENTS 



We estimate Bogoliubov coefficients in order to define Hawking temperature for a nonstationary background in this 
section. We consider an eternal black hole in thermal equilibrium with its surroundings, that is, in the Hartle-Hawking 
state. In asymptotically flat cases, we should immerse the black hole into a thermal bath to realize this state. In 
asymptotically AdS cases, on the other hand, it is naturally realized due to its boundary condition [2^, because 
the black hole is enclosed in the AdS boundary. 

We exploit the saddle point approximation to evaluate the Bogoliubov coefficients and determine the temperature 
from them. After reviewing usage of the approximation for the static spacetime in Sec. Ill A[ we consider its extension 
to nonstationary spacetime in Sec. IIIBl In this extension to nonstationary spacetime, we introduce a quantity, k{u), 
which may be interpreted as an extension of the surface gravity of static Killing horizon. In Sec. Ill CI we reinterpret 
this quantity from a geometric point of view, and clarify that k(u) is naturally associated with the past horizon of 
the eternal black hole. 



A. Static spacetime and the saddle point approximation 

For simplicity, we consider the two-dimensional part of the spacetime consisting of time and radial directions, 
that is, we focus on the s-wave sector of radiation. We also adopt the geometric optics approximation, and neglect 
the backscattering of the waves due to the curvature of spacetime. We introduce the null coordinate u which gives a 
natural time for observers in the asymptotic region (null infinity). If the spacetime is stationary, this time corresponds 
to the Killing time. We introduce another null coordinate U which is the affine parameter on the past horizon. For 
the stationary case it becomes the familiar Kruskal coordinate. Because lines described by u = const, are outgoing 
null geodesies, geodesic equations give a relation U = U{u) between the two null coordinates. 

Now, we consider a massless scalar field. In the current setup, solutions of the field equation are simply given 
by arbitrary functions of each null coordinate. In a standard manner we can define positive frequency modes with 
respect to u and U , respectively. Then, the Bogoliubov transformation between two sets of modes cx e~*"" and 
ua, cx e"'"*^ is determined by the Bogoliubov coefficients 

""4 = ±A,/^ r d^— e±-«-*"^("), (1) 
PujCj J 27r V w J_oo du 



^ The essential part of the calculation is similar to Ref . Q . 

^ The derivations are quite similar even on the background spacetime with gravitational collapse and black hole formation. Sec also 
comments in Sec. IIVI 

^ In asymptotically AdS cases we should consider a conformally coupled scalar field. 
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where the upper and lower signs correspond to a^^j and /3^un respectively. 

Now, to evaluate those coefficients by the saddle point approximation, we will consider the integral 

exp(/)(u)du, (2) 

^ — oo 

where 

(/)(it) = log — — ± icju — iCjUiu). (3) 
au 

Saddle points are located aX u — u^, given by = 0, where 

(j) [u) = —log — ±iuj - luj—, (4) 
au au au 

and a prime denotes a derivative. 

Consider a static black hole, for example. When the surface gravity at the Killing horizon is k, the relation between 
the two coordinates u and U is given by U{u) = — exp(— km), which is the coordinate transformation for the Kruskal 
coordinate. Then, we have 

= -K ± - iwKe"''"' = 0, (5) 

and the saddle point 

= -K-ilog— =— log-r-- -T^ , 6 

where re^** = 1 ^ iut/n. We notice that the real part of depends on u} and it becomes larger as uj increases. 
This implies that the Kruskal modes with very high frequencies lo are relevant for late-time features, while relevant 
frequencies uj will change depending on the time of observation. [See also discussions around Eq. pO[) .] Using 



(t){u^) = log K + re^^" log 



(7) 



we have the saddle point approximation of the integral ^ as 



("Mm e^(-) / e^(— )Mm = e^(-\/-^ = « { ^^^"l (8) 



which becomes a good approximation when lo/ 1. This expression derives the Bogoliubov coefficients of Eq. ([T]) 
given by 

exp . (9) 



I ^UJU) p V ^ 

Thus, we have the familiar result of the Hawking temperature to be T = k/2'k by the saddle point approximation. 
We could also perform the integral explicitly to obtain a well-known exact expression [l[ 

e^(")dM = K(-^w)-'■^^"'^ (re^^^) . (10) 

Using the Stirling formula of the gamma function 

r(i + z) = zr(z)~y2^(^)\ (11) 



we have Eq. ([8]) again from Eq. ([T0| . This result indicates that when the geometric optics approximation is satisfied, 
namely, when uj/ 1, the saddle point approximation is valid. 
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B. Extension to nonstationary spacetime 

Let us now take a nonstationary spacetime background. To probe this spacetime, we use wave packets which 
are locahzed in both the time and frequency domains rather than plane waves spreading over the whole time. The 
wave packet is peaked around a time u = Uq with width ~ Au, which are the time and duration of the observation, 
respectively. The Bogoliubov coefficients for a wave packet arc obtained by inserting a window function into the 
integrand of Eq. ([T|), that is, 

where wau{x) is a window function which goes sufficiently fast to zero outside the interval Au around x = 0. In 
addition, we assume that analytic continuation of WAuix) is varying slowly at least for |x| < Aw in the complex plane. 
Using the Klein-Gordon product, these coefficients are written as A^j^} — {iPcj,uq) and Bi^^ — —(ip^,Uuj), where we 
have constructed the wave packet localized around time u = uq and frequency w as 

^^(u) = -^WAu{u - uo)e-*"". (13) 

We note that A^cj and B^^cj satisfy the following relation 

da; (l^^^p - ^ (V>„,V'^), (14) 

which follows from the completeness of {ucj}. The number density in terms of the wave packet mode ipi^ is given by 

1 



(V'a;,V'c^) Jo 



d^lB^^^r, (15) 



which is roughly the number of particles generated in the frequency band of width ~ 1/Au around frequency w. 

Let us suppose that the evolution of the background spacetime is sufficiently slow around u ~ uq within time 
interval Au. To clarify what slowly evolving is, we shall introduce the following quantity 

/ N d , dC/ 
,(,)^__log_. (16) 

If the background spacetime is stationary, k(u) becomes a constant and it is nothing but the surface gravity of the 
black hole (Killing) horizon, as we have seen in the previous section. Therefore, we may state that the spacetime is 
evolving slowly when k(u) is almost constant during Au. In that case, we may express k(u) as 

k(u) = Ko[l + e/(u)], (17) 

where /(u) is a real function which satisfies /(uo) = and |/(m)| < 1 for |u — uq\ < Au, and e is a small parameter. 
In other words, this assumption implies 

An 

<e for |u-uo|<Au, (18) 

Ko 

where Ak = |k(u) — ko|. For discussions below, we further assume that analytic continuation of /(u) satisfies |/(u)| < 1 
for any complex u such that |u — uo| < Au. 

Then, we may expand U'{u), U{u) and U"{u) with respect to e as 

U'{u) = Uq exp ^- J K{x)dx^ = Uq exp {-kq {Su + eg{u))} = UqC''^"^'' (l - enog^u) + 0{e^)) , 

U{u) = C/o + / U'{x)dx = C/o + C/q EKo / e-««("-"'')g(a;)dx + 0{e^) , 

U"{u) = -k{u)U'{u) = -Ko (1 + e/(u)) C/^e-'^«^" (l - €Kog{u) + 0{e^)) , 
where we defined C/q = U{uf)), Su = u — uq and 

5(u) = f f{x)dx. (20) 



(19) 
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Note that \g{u)\ < \6u\ follows from the assumption < 1 for \Su\ < Au. We also expand the saddle point u = 

as 



,(0) 



,(1) 



(21) 



Then, the equation for the saddle point from Eq. (jU is expanded using Eq. (fT9)) as 

= = -Kore"^'' - iwC/^e-««''"-'" + cko [-/(ul"') + ituU^e-''"'^'" + g{uf^))] + O(e'), (22) 

where we have redefined re^'^ = 1 =p iw/kqi which satisfies r 3> 1 because of the geometric optics approximation 
lo/kq ^ 1. Solving 0'(w,) = order by order of e, we find 



6u 



(0) 



Kn ^ loS 



iKQTe 



= — — log 



,(1) 



(23) 



Kore 



Plugging Eq. ([23| into Eq. ([3]) and (/'"(u) = —k'{u) — iCjU"{u) at u = u, and expanding with respect to e, we find 
0(u*) — 0(wo) = '/'o + e*/*! and (^"(w*) = <^g + e<^" are given by (/)(mo) = log ± iwug — «wJ7o, 



«)o = -re 



(ko'5"1"' + l) 



g{u)du 



and 



-Kr,re 



(24) 



(25) 



where ©(e^) terms are omitted. 

To guarantee that the perturbative expansion above is valid and that the calculation of the saddle point approxima- 
tion is not affected by the correction terms, we should require each of correction terms [eui^\ e4ii and t4>'l in Eqs. 



(HH) and ([^5]) ] to be much smaller than their leading term (5u'^\ and 0q, respectively). As for e|u*'^''| <C \5u 

and e|(/>i| ^ l^ol; can see from Eqs. (|23p and ([^5)) that we should assume both e ^ 1 and <C kq?" to 

hold to satisfy such requirements, where we used \g{u)\ < \Su\. As for e|Re(/)i| ^ |Re(/)o|, we have 



,(i)| 



(0)1 



Re 00 = T wlmuj, ' — KoKeSui — 1, 

(/>Rc u 
Jrc «1' 



Rcu'-">+ilmu'-" 



-'^-'^g{u)du , 



and 



Re(/)i| < Kolf^ul"''! + KQr|(5M*''' I llmu!;' 



(0)1 



,(o)| 



(26) 



(27) 



Then, for e|Rc0i| ^ |Rec/)o| to hold, we need to additionally require eKolf^itl^' | <C 1. Under these assumptions, we 
find that the integral in Eq. (jl2p can be evaluated by the saddle point approximation as 



,(0)^ 



2n 



\<f>l' 



(28) 



So that the suppression due to the window function, w^u, is weak and the integral value becomes close to that for 
the unwindowed function, the saddle point should satisfy |(5m,| < Am. For the imaginary part of itl"'' the condition 
|Ju»| < Am becomes 



1 
Au 



< 



^0, 



and for the real part it becomes a condition on cj, given by 



Kor 



(29) 



(30) 
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If Cj is out of this region, no saddle point exists between |u — uqI ^ Au. and then the integral will be suppressed due to 
the window function. This condition means that only the Kruskal modes of a) with limited frequency band, specified 
by Eq. (|30p . can have correlations with the wave packet mode ipuj due to its localization in the time domain. 
Consequently, sufficient conditions for the saddle point approximation to be valid arc given by 

cKo < -r^- < Ko < w, \f'{u)\<^KQr, (31) 

for |u — wqI < Am. Note that the condition on e, ekq ^ Am~^, can be rewritten as AkAu <|; 1, where Ak is defined 
by Eq. ([T5| . If we can take the time interval Au(> ^) satisfying this condition, we have a sufficiently small e ^ 1 
for which Eq. ([5T|) holds. Roughly speaking, this condition is indicating that the time variation of k{u) should be 
sufficiently moderate to satisfy ^ ^ between the interval |m — uo| ^ . Under these assumptions, at the 

leading order we have 

which implies the spectrum observed at the time u = uq becomes a thermal one with temperature Ko/2n for high 
frequencies w 3> kq. 



C. Surface gravity for past horizon 

In this section, we consider the geometrical meaning of k{u) defined by Eq. (|16p in the previous section. Ingoing 
null vectors with respect to null coordinates u and U arc written as 

dV r„ fd 



respectively. Note that U is an affinc parameter on the past horizon because it is one of the Kruskal coordinates. 
Using fc" = U'{u)k°', we have 

fc'^Vafc'' = rVjL/'(M)P) = (fc"VaJ7'(u))fc'' + {U'fk''Vak\ (34) 

The last term vanishes because U is the affine parameter. As a result, we obtain k°-V ak^ = —n{u)k^, where 

'^(w) = -T"log^- (35) 
au au 

We shall call k{u) "surface gravity" for the past horizon because it describes the inafhnity of the null generator 
of the past horizon, which is defined by the asymptotic time at the null infinity. When the spacetime is stationary, 
the past horizon and k{u) coincide with the Killing horizon and its surface gravity, respectively. In this sense, niu) 
is a natural extension of surface gravity of stationary spacetimes. (In Ref. Q, it was called "peeling properties" for 
gravitational collapse.) It is worth noting that the surface gravity for the past horizon is not determined only by 
local geometrical quantities on that horizon but also by the relation between the time coordinates on the horizon and 
in the asymptotic region. In fact, the Killing vector which defines the surface gravity of the Killing horizon is also 
determined using the asymptotic time in a stationary case. 

Roughly speaking, the particles of the Hawking radiation observed in the asymptotic region start from the vicinity 
of the horizon and propagate along outgoing null geodesies. Every outgoing null ray arriving at the null infinity expe- 
riences the near-horizon region of the past horizon rather than the (future) event horizon. Based on this observation, 
it seems to be natural that the spectrum is affected by the surface gravity near the past horizon. 



III. APPLICATIONS 



In this section, we consider explicit examples of transitions from an initial equilibrium state to another equilibrium 
state. First, as a simplest example, we consider null shell accretion into a black hole in asymptotically flat spacetime. 
Next, we will focus attention on the AdS-Vaidya spacetime. This is a toy model to describe a thermalization process 
in the AdS/CFT correspondence. 
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A. Asymptotically flat case: null shell accretion 

111 this section, we consider a black hole in asymptotically flat spacetime, whose mass is initially Mi and changes 
into the final mass Mf due to accretion of a null shell. Then, the Hawking temperature should initially be that of 
the initial black hole and eventually become the temperature of the final one. The metric of a static and spherically 
symmetric black hole solution is given by 

dr^ 

ds2 ^ ~fi{r)df + + r^dn\ (36) 

where the subscript / is i or /, which are for initial and final quantities, respectively. Double-null coordinates (u, v) 
are given by 

dr 

u + v = 2t (37) 



fir)' 

To describe null shell accretion, we will join two spacetimes at = in a manner such that for v < we use the 
metric of the initial black hole and for v > that of the final one. At v = the junction condition implies that the 
radial coordinate should be identified as 

/dr f dr , , 

77T, "f = -2/7TT- (38) 
Mr) J fi{r) 

Also, it leads to the following relations 

dR^_MR)_ di?^_ifCR) .^q^ 

dUi 2 ' duf 2 ' ^ ^ 

where R is the radius of the null shell, and these are nothing but the equations of motion for the null shell. 
The Kruskal coordinate at the past horizon, namely, that of the initial black hole, is 

[/ = - exp(-KiUi), (40) 

where Ki denotes the surface gravity of the initial black hole, and the retarded time at null infinity is that of the final 
black hole as u = U{. 

dU , , dui , , /f (i?) , ^ , 

— = Kiexp(-KiUi)— = Kiexp(-KiUi)-— — (41) 
du du Ji(R) 

The surface gravity of the past horizon is 

Kin) ^^—\og—^K-^- IMl^ + //(^) dR 
du du ' du fi{R) du fi{R) du 

mR)\ ff{R) , m) ^ ' 



MR) 



where R{u) denotes the radius of the null shell at time u. For early time u — ?> — oo, the radius of the shell becomes 
i? — >■ oo and we have k{u) — ^ Ki. For late time u — > oo, the radius of the shell approaches the horizon radius of the 
final black hole R — > rt and we have k(m) — > Kf as u — > oo, where we have used Kf = ff(r{)/2. These imply that the 
asymptotic observers detect a change of the Hawking radiation at the retarded time when the null shell comes into 
the vicinity of the black hole horizon. 

Consider the four-dimcnsional Schwarzschild case with /i(7') = 1 — 2Mi/r for example. We have 

_ R{u) - 2AM d AMjRju) - 2M,) 

^^''>- iMiRiu) ' du^^""'- AMiR{u)^ ' ^^^^ 



where AAf = Mf - Mj. If Mi > AM, the condition 

1 dK ^ AM 



k2 du ^ Mi 



< 1 (44) 



is satisfied for any retarded time u. Therefore, it is concluded that the Hawking temperature observed by asymptotic 
observers at time u is given by k(u)/2tt in the current case. We note that k{u) changes gradually even though the 
spacetime describing the null shell accretion is not smooth. 
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B. Asymptotically anti-de Sitter case: AdS— Vaidya 



Next, we discuss an asymptotically anti-de Sitter background. In this case, null infinity at which the Hawking 
temperature would be observed is not a null surface but a timclikc surface, namely, the so-called AdS boundary. 
Moreover, even when there is no black hole horizon, asymptotically AdS spacetime has a past horizon (not a white 
hole horizon but a Cauchy horizon). 

We consider the 5D AdS-Vaidya spacetime 

ds^ = ^[-F{v, z)dv^ - 2dvdz + dxg], (45) 

where F{v, z) is given by 

F{v,z) = I -2m{v)z^, (46) 
and the curvature radius of the AdS is set to unity. Wc suppose that the mass function m(v) is 

mo {v < 0) 

m{v) = ^ jno + Amsin'^ 2A~ (0 < < Av) . (47) 
mo 4- Am [v > Av) 

During an interval Aw, the null fluid is injected into the black hole with mass mo, and the mass eventually becomes 
mo -|- Am. For the AdS-Vaidya, Bondi energy observed at the boundary is described by the mass function m{v). 
Note that the time coordinate v is nothing but asymptotic time at the boundary. In the context of the AdS/CFT 
correspondence, m(v) corresponds to the energy density of the CFT matter. Therefore, the duration of the injection 
Av will represent the time scale of energy density change. 

Now we shall consider the Hawking temperature in the current case. We introduce double-null coordinates (w, v) 

as 

z — z{u,v), V ^ V. (48) 

In these coordinates, the boundary 2 = is described hy u = v. An asymptotic time t at the boundary is given by 
t = iii^, and it is equal to (or u) there. The coordinate condition for {u,v) to be the double-null coordinates leads 
to 

dz _ F{v,z) 
dv~ 2 ' 

which is equivalent to the geodesic equation for the outgoing null ray described by u-constant line. Note that we have 
^1 „ = — 4^1 „ = — ias the boundary condition at z = 0. 

dv 1 2=0 au lz=0 2 ^ 

In order to calculate the surface gravity for the past horizon we need to know the relation between the asymptotic 
time and the affine parameter at the past horizon. Because the spacetime is static (strictly speaking, independent of 
v) for the initial black hole region (v < 0), the canonical null coordinate u is given by 

2 

du = dv + —dz, (50) 

where we note that the metric function F depends only on z. Moreover, the afhne parameter at the past horizon, 
namely, the Kruskal coordinate U is immediately given by U(u) = — cxp(— Kju), where Kj is the surface gravity 
determined by the initial black hole with the mass mo. It is worth noting that, in general, the above null coordinate 
u is different from the asymptotic time u defined previously. Let us recall the definition of u. The differential form 
becomes 

dz dz F 

dz = -T—du + —dv — (— Adu + dv), (51) 

ou av 2 



where A is defined by 



2 dz du 
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Boundary time (i(=v) Boundary time (m=v) 



FIG. 1: K,{u) as a function of time {u = v) at the boundary (left) and behaviors of k{u) near the final temperature Kf (right). 
Time scales of injection are taken to be Aw = 0.1, 0.2, 0.5, 1, 2, which are shorter than a typical time scale given by temperature. 
We set mo — 1/2 and Am — 15/2, which give ki = 2 and Kf = 4, respectively. 



If the spacetime is static, we can take the canonical double-null form such that A = 1 everywhere. It turns out that 
A describes the relation between the asymptotic time u at the boundary and the canonical time u in the initial black 
hole region, that is, the redshift factor for the outgoing null ray which goes from the initial black hole region to the 
boundary, due to dynamical background. (See the Appendix for details.) 

Now, we can obtain A(m) as follows. Differentiating Eq. with respect to u, we have 

d dz 1 dF dz 
dv du 2 dz du 

Integrating the above equation and the geodesic equation with respect to v from the boundary v = u io the initial 
static region {v < 0), we obtain |^ for w < 0. Note that we must integrate those to the past horizon {v = — oo) in 
general. However, in the current case it is enough to evaluate only up to w = 0, because the spacetime becomes static 
for < and the redshift factor A(u) does not change any more. 
As a result, we have 

k{u) = log ^ = KiA(u) - log A(u), (54) 
du du au 

where we have used 4^ = A4^. 

du du 

In Fig. [T] we plot k{u) as a function of the boundary time [u = v) for shorter injection times Av = 0.1, 0.2, 0.5, 1, 2. 
The initial and final temperatures are given by Ki = 2 and K{ ~ 4, respectively. The conditions of Eq. (|3T|) are almost 
satisfied for these parameters and thus the results of the saddle point approximation are valid. ^ In these cases, we 
find that k(u) converges into the final temperature exponentially after a transient phase. Time scales of the whole 
process of temperature change are roughly given by ^ 3 and do not strongly depend on time scales of injection 
but are dominated by exponential relaxation to the final temperature. Recalling that the time scale determined by 
temperature T of a black hole is given by 1/T — 27r/K, we may interpret this result as the time scale of the temperature 
change is governed by the temperature of the final black hole. See also Fig. [2] about final-temperature dependence of 
relaxation rate for the case of null shell accretion, for which the injection is instantaneous and equivalent to the case 
of Av = 0. Note that since we are focusing on time scales shorter than the thermal time scale Am ^ k^^ and the 
minimum time resolution of observation is given by this thermal time scale, we should interpret k(u) not as a value 
at time u but rather a time-averaged value over time scale ^ k~^. 

As the injection time becomes longer than the thermal time scale, k(u) begins to show qualitatively new behaviors. 
In Fig. [HI for Aw = 0, 2, 10, 50 we plot k{u) together with a quasistatic temperature 27rrqs = 2 (2m{v))^^^ , which is 



If the difEcrcncc between the initial and final temperatures becomes very large, the adiabatic condition may be temporarily violated. 
This means that the system is so far from equilibrium that temperature can not be defined at that time. However, after the system has 
relaxed near the final temperature, the interpretation of k{u) as temperature becomes well-defined again. 
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decay phase is proportional to inverse of the final temperature. 
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FIG. 3: k{u) and a quasistatic temperature 2nTcis = 2(2m(w)) ' for various injection times Av 



naively determined from the mass functions m{v) in a similar manner to static black holes. We can observe that both 
curves begin to coincide as the injection time becomes longer. In other words, the variation of the temperature k{u) 
tends to follow the variation of the black hole mass m(v). 

These results lead us to the conclusion that, even if the mass injection is rapid, the temperature will gradually 
change over a time scale no shorter than the thermal time scale ^ In this sense, we may say that there is a 

finite relaxation time irrespective of the speed of energy injection. When energy injection and the resultant black hole 
evolution is slower than this relaxation time, the temperature tends to evolve in a similar way to the mass injection. 
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IV. SUMMARY AND DISCUSSION 

111 this paper we have discussed Hawking temperature for nonstationary spacetimes. We introduced a useful measure 
for time evolution of spacetime, and showed that the Hawking temperature at each time is determined by the surface 
gravity of the past horizon when the time evolution of the spacetime is sufficiently slow. This definition of temperature 
is quite natural because observers can determine it without using causally inaccessible information, such as the position 
of the event horizon. 

Although we have considered an eternal black hole which has the past horizon, the essential calculation does not 
change for black hole formation in an asymptotically flat spacetime. In this case we should draw null rays from 
future null infinity back to the past null infinity instead of the past horizon, and as a result we can reproduce the 
original derivation of the Hawking radiation. The definition of k{u) is the same as Eq. ()16p . even though it cannot be 
interpreted as "surface gravity" and should be understood rather in the context of "peeling properties" of outgoing 
null geodesies, which was discussed in Q. It would be interesting to study the properties of k{u) for black holes 
in an expanding universe, such as those discussed in Refs. [sl, [2^, for which the derivation is similar to the case of 
asymptotically flat spacetime. In the case of asymptotically AdS spacetime, in the Poincare chart the past horizon, 
which is the Cauchy horizon, exists even in the pure AdS spacetime. This setup is quite similar to the eternal black 
hole, and we can naturally interpret k{u) as the "surface gravity" of this past horizon. 

In the context of the AdS/CFT correspondence, our results indicate that quantum fluctuation of the boundary CFT 
is governed by the surface gravity of the past horizon when the bulk is the eternal black hole, namely, the thermal 
equilibrium state. Particularly, our result can be applied to phenomena discussed in Refs. [l^, Ell, which are directly 
related to the Hawking radiation in the bulk. It may be interesting to study its implications further. It would be 
also interesting to study k{u) in numerical solutions obtained in, e.g., Refs. (l3. [l5l. [23j . though in general numerical 
calculation of k{u) becomes difficult at late time due to exponential pileup of outgoing null geodesies onto the horizon. 

It is shown that the change of the Hawking temperature is delayed by thermal time scale compared with the mass 
injection from the boundary. This behavior can be interpreted as follows: The injected matter falls down into the 
black hole and changes the geometry in the vicinity of the horizon. After that, the Hawking radiation affected by 
this change emanates from there and reaches the boundary. This process of infall of the matter and return of the 
modified radiation will take finite time and, roughly speaking, this is the origin of the delay we observed. From the 
viewpoint of the boundary theory, this delay time can be interpreted as the relaxation time scale needed to achieve 
thermal equilibrium after the energy injection. 

For planar AdS black holes (namely, large black holes in asymptotically AdS spacetime) the horizon is located 
relatively near to the AdS boundary. If a black hole evolves sufficiently slower than the delay time, we may neglect 
its delay. In that case, the quasistatic approximation, such that one synchronizes dynamics of the event or apparent 
horizon and the boundary in terms of the advanced time and uses quantities associated with the future horizon, might 
be justified. 

What we emphasize is that the "surface gravity" of the past horizon, is not only a geometrical quantity but 

also one with a physical meaning in the sense that it governs the thermal spectrum of the Hawking radiation observed 
by asymptotic observers. Many of thermodynamic properties of black holes, however, tend to be associated with the 
future event or apparent horizon in the previous works, not with the past horizon. It is open to discussion how these 
points of view are related to each other. It may be also fruitful to study how k{u) is related to other probes of the 
black hole spacetime, such as those discussed in Ref. [l^, and how our approach is related to other derivations of the 
Hawking radiation, such as the tunneling approach. 
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Appendix A: redshift factor 

111 this appendix we show that X{u) defined by Eq. (|52p is the redshift factor for the outgoing nuU ray. Now, we 
recaU that the two-dimensional part of the metric is given by 

ds^ = ^\-F(v, z)dv^ - 2dvdzl (Al) 

A coordinate transformation such as z = z(u, v) and v = v to the double-null coordinates {u, v) gives us the differential 
form 

dz d z F 
dz= —du+ —dv ^ ~ — (-Xdu + dv), (A2) 
au ov 2 

where we have defined as A = — ^ / 

OU I ov 

We consider a timelike vector field = d/dv which characterizes a natural Killing time in stationary regions and 
also in asymptotic regions near the boundary. Note that frequencies of the Killing modes are defined with respect to 
this Killing time. By using the 1-form it is rewritten as 

^a = -^iFdv + dz). (A3) 
z^ 

The tangent vector of outgoing null geodesies is given by 

r EE A, (A4) 
ds 

where s is an affine parameter, and then the geodesic equations lead to 

d t, OF f . F ^ 

ds z'' az 2z^ 2 

where the dot denotes the derivative with respect to s. We note that the last equation is equivalent to the coordinate 
condition dz/dv = — F/2 defining the double-null coordinates. 

The "Killing energy" E associated with the null geodesies is defined by 

E^-na^^{Fi + z) = ^^. (A6) 

If is truly a Killing field, E should be constant along the null geodesies. From the geodesic equations we have 

d V 1 dF V ^^^^ 
dv z^ 2 dz 

where we have used d/ds = vdjdv. On the other hand, differentiating dz/dv = —F/2 with respect to u, we have 

d dz 1 dF dz 
dv du 2 dz du' 



As a result, we find 



'/ = ci, (A9) 
du V 

where C is an integration constant which can be absorbed into normalization of the affine parameter. Then, the 
"Killing energy" can be rewritten as 

dz / dz 

dv / du ^ ^ 



Since the boundary conditions at z = are given by 
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wc have -Eb = which is the energy observed at the boundary, and 

EU^^is) = ^b/A. (A12) 

Using the geodesic equations, we also have 

dE _FdFv' V dF _ V ^dF _ .d\nF 

di" ~ 2" 5^2^ 2^ d7 ~ 2^^9^ ~ ^ ' 

Integrating it along the outgoing null geodesies described by 2 = Zg{v), we obtain 



E^=E,e^vl dv[^^] , (A14) 



dv 



where Ei and -Eb are the Killing energy observed at an initial surface and the boundary, respectively. 

Consequently, for the outgoing null ray described by u = const., we have the redshift factor from an initial time Vi 
to asymptotic time u = v^, as 

X(u)=e^p r dv(^^] . (A15) 

It turns out that if the spacetime is stationary, namely, F does not depend on v, we have A = 1. The redshift factor 
X{u) deviates from the unity when time dependence is turned on. 
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